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Abstract 

We study the possible origin of hydrodynamic turbulence in cold accretion disks 
such as those in star-forming systems and quiescent cataclysmic variables. As these 
systems are expected to have neutral gas, the turbulent viscosity is likely to be hy- 
drodynamic in origin, not magnetohydrodynamic. Therefore MRI will be sluggish or 
even absent in such disks. Although there are no exponentially growing eigenmodes 
in a hydrodynamic disk, because of the non-normal nature of the eigenmodes, a 
large transient growth in the energy is still possible, which may enable the system 
to switch to a turbulent state. For a Keplerian disk, we estimate that the energy 
will grow by a factor of 1000 for a Reynolds number close to a million. 
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1 INTRODUCTION 



In the context of Keplerian disks, while hydrodynamic simulations seem to 
confirm the absence of turbulence in the non-linear limit [1], the Magneto- 
Rotational Instability [2] , discovered by Balbus and Hawley [3] within ionized 
accretion flows, confirms the existence of MHD turbulence. This helps to un- 
derstand the origin of shear stress or "viscosity". However the MRI driven 
instability fails to operate in disks with small ionization fraction. Examples of 
such systems are proto-planetary disks, outskirts of AGN accretion disks [4] 
etc.. As a result, the route to turbulence and subsequently accretion in neutral 
disks has remained one of the outstanding puzzles in modern astrophysics. 

Laboratory experiments of Taylor- Couette systems (which are similar to ac- 
cretion systems) seem to indicate that, the flow is unstable to turbulence for 
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Reynolds numbers larger than a few thousand [5], even for subcritical sys- 
tems. Based on this, Longaretti [6] concludes that the absence of turbulence 
in previous simulations [1] is due to their small effective Reynolds number. 

The fundamental reason for subcritical turbulence may be the result of definite 
frequency modes not being orthogonal in a shear flow. Therefore, a suitably 
tuned linear combination of modes can show an arbitrarily large transient 
energy growth 1 . This transient growth may possibly lead to sustained hydro- 
dynamic turbulence for large enough Reynolds numbers. 

Here we analyze this transient growth phenomenon in the framework of ac- 
cretion disks. We present the set of scaling relations for the transient growth 
as a function of the Reynolds number. 



2 BASIC EQUATIONS AND SCALING RELATIONS 

Let us define the shear frequency 2A and the vorticity frequency 2B as, 2A = 
—qfl, 2B = (2 — q)Q. Then the momentum balance equations give 

£ = Mto - f + ,W % = -2Bu - f + W\ d J = -f + *W (1) 
at ox at ay at oz 

along with the incompressibility relation du/dx + dv/dy + dw/dz = 0, where 
p = P/p. Symbols u,v,w and Q are three velocity components and angular 
frequency respectively, v, P, p are dynamical viscosity, pressure and density of 
the system respectively and q is a parameter which controls the flow pattern 
(whether a Keplerian disk, a constant angular momentum disk or plane Cou- 
ette flow). The Lagrangian time derivative d/dt is, d/dt = d/dt — qQxd/dy. 
Here x varies from — L to +L, and y and z from to 2rc/k y and 27i/k z respec- 
tively, where k x , k y , k z , are the corresponding components of wave- vectors. 
The boundary conditions are u = v = w = 0aix = ±L. 

2. 1 Constant Angular Momentum Disks and Plane Couette Flow 

It is clear from (1) that the set of equations is very symmetrical for a constant 
angular momentum (q = 2) and plane Couette (Q = 0) flow, except that x 
and y are interchanged. Therefore the growth is identical in the two cases. As 
the previous analytical/numerical studies (e.g. [8,7]) showed that the growth 

1 Growth is defined as the ratio of energy at a particular instant to that at the 
beginning. 
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is maximum for vertical perturbations (k y = 0), here we consider such pertur- 
bation. Noting that the Reynolds number is, R = \1A\L 2 jv, we find that the 
maximum growth in energy for a given k z is 



' u 2 {t) + v 2 (t) + w 2 (t) \ 



R 2 k 2 L 2 e- 2 



±vr 2 + k 2 L 2 



Maximizing this over k z , we find that the optimum wavevector is k z L 
1.57. The maximum growth factor and the corresponding time are 



(2) 
tt/2 = 



Gmax (R) = 0.82 x 10- 3 R 2 , \2A\t max {R) = 0.13/2. 

(3) 

Clearly the growth scales as R 2 and therefore can be very large even for modest 
R. For example, when R = 1200, the above relations give G max (R) > 1000, 
which is probably large enough to induce non-linear feedback and turbulence. 



2.2 Keplerian Disk 



Here we study eqns. (1) for q = 1.5. As the earlier analysis [7] hinted that 
k z = generates the best growth, we consider such perturbations. 

We consider a plane wave that is frozen into the fluid and is sheared along with 
the background flow. If the flow starts at time t — with initial wave-vector 
{k X i, k y }, the k x at later time is given by k x (t) = k X i + qVtk y t. The maximum 
energy growth is then approximately 

G™A„ *„ R) ~ (17) ff^ exp ■ (4) 

when k x>min L = 1.7. Maximizing this with respect to k xi and k y , we obtain 

G m ax(R) ~ 0.13i? 2/3 , |2A|t ma;r (i?) ~ 0.88R 1/3 . (5) 
Also maximizing with respect to k xi , keeping k y and R fixed, we obtain 



G ma x(k y , R) — 27 — pr^-e 2//3 i? 2 ^ 3 , |2A|t maa; (A;j / , i?) — (k y L) 2 ^R 1 ^. (6) 

We see that the maximum growth scales as R 2 ^ 3 . Though this rate of increase 
is less than that for a q = 2 disk and plane Couette flow, still, at a large enough 



3 



R, the growth can become large, and may cause turbulence in a Keplerian disk. 
For example, for R = 10 6 , G max (R) > 1000, which by analogy with a q = 2 
disk may be enough to cause turbulence. 



3 CONCLUSION 

We have shown that significant transient growth of perturbations is possible in 
a Keplerian flow. Although the system does not have any unstable eigenmodes, 
because of the non-normal nature of the eigenmodes a significant level of 
transient energy growth is possible at a large Reynolds number for appropriate 
choice of initial conditions. We argue that a plausible critical Reynolds number 
for sustaining hydrodynamic turbulence in a Keplerian disk is ~ 10 4 — 10 6 (for 
detailed discussions, see [7]). We base this on the expectation that once the 
growth crosses the threshold, secondary instabilities of various kinds, such as 
the elliptical instability (see e.g. [9]), might serve as a possible route to self- 
sustained turbulence. However, it remains to be verified that these instabilities 
are present. Also, even if they are present, one will need to investigate whether 
they lead to non-linear feedback and 3-dimensional turbulence. 

This work was supported in part by NASA grant NAG5-10780 and NSF grant 
AST 0307433. 
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